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Abstract 

I give a short review of the relation of infrared renormalons in QCD and higher twist effects, with 
the emphasis on possible applications. In particular, I present estimates of renormalon-induced 
uncertainties in deep inelastic sum rules and explain how the renormalons can potentially be used 
to unravel the structure of nonperturbative effects in complicated situations and to indicate possible 
systematic sources of large perturbative corrections. 
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1. The modern precise data on "hard" processes in QCD require a theoretical descrip- 
tion to power-like accuracy. A clear example is provided by the Gross-Llewellyn Smith sum 
rule in the deep inelastic scattering, where the higher-twist (HT) 1/Q 2 correction produces a 
major uncertainity in determination of a s [|lj. For a generic physical observable dominated by 
short distances one expects a theoretical prediction of the form, schematically, 

^)=^ e [l + r 1 ^ + r 2 (^) 2 + ...] + E^, (1) 

where a s = a s (Q) and H s 's are dimensionful nonperturbative parameters (with dimension 
2s > 0) describing the HT corrections. A conceptual problem, which I am going to review in 
this talk, is that the discrimination between perturbative corrections and HT contributions is 
ambiguous. This will imply that it is not possible to attribute a fully quantitative meaning to 
power-suppressed corrections, unless some prescription is used to sum the perturbative series. 
On the other hand, ambiguities in summation of the perturbative series can serve to indicate 
which powers of 1/Q are required in the sum in ([!]). 

Assuming for simplicity a (euclidian) quantity dominated by a large scale Q, the leading- 
order correction involves the gluon exchange with the large virtuality of order Q. Progressing 
to a higher order n the average gluon virtuality still remains proportional to Q, k ~ a n Q where 
a n is a certain coefficient, simply because there are no more dimensionful parameters. However, 
the a n can (and do) decrease with n, so that in very high orders n such that a n ~ Q/Aq CD 
the perturbative calculation fails. An inspection shows that the most dangerous Feynman 
diagrams are those related to running of the QCD coupling, in which the average virtuality 
decreases exponentially a n ~ exp[— n/s], where s is some number. It is possible to show 
that this intervention of infrared (IR) regions reveals itself in a rapid - factorial - increase of 
perturbative coefficients in high orders ||: 

R(Q) = i^ ee £r n <(Q); r n ~ ((3 /s) n n\ (2) 

n 

where (3q = (11 — 2/3n/)/(47r). The physical origin of the large coefficients is simple: the 
gluon exchange with virtuality k involves the QCD coupling at this scale a s (k). However, the 
perturbative expansion (fj) is assumed to be in powers of a s (Q) at the scale of the external 
momenta. Thus, we get large coefficients simply by reexpressing a s (k ~ exp(— n/s)Q) in 
terms of ot s {Q). 

A factorial growth of perturbative coefficients means that the perturbative series is at best 
an asymptotic series: the fixed order contributions r n a™ decrease at small n, reach a certain 
minimum value at n = n ~ l/a s , but then again start to increase and blow up. This means 
that a perturbative calculation only makes sense up to the order n Q , and the accuracy of this 
calculation, or, equivalently, the effect of the "tail" with n > no is of the order of the minimum 
term 

oo \s 

£ r n a n s {Q) ~ r no aT(Q) ~ exp[-s/(/WQ))] ~ "Sr ' ( 3 ) 

n=n \ V / 

where I used the asymptotic form of the coefficients (H) and the one-loop formula a s (Q) = 
1/(A) m (Q 2 M 2 ))- 111 the common jargon, the divergences of perturbative expansions are called 
"renormalons" [|J, and I will refer to the ambiguity in the summation of the series (0) as to 
the "renormalon ambiguity". 

The deficiency of the perturbation theory has a profound reason, indicating that calcula- 
tion of physical quantities to power-like accuracy requires taking into account nonperturbative 



effects. The renormalon ambiguities must be compensated by ambiguities in the HT correc- 
tions complementing truncated perturbative expansions. Hence, the required values of s in 
(§) have to be in one-to-one correspondence to the required HT corrections in ([!]). This im- 
plies that the required powers of 1/Q can in principle be determined from purely perturbative 
calculations. In practice, already a simple approximation (referring to the 1/Nf expansion) 
usually recognizes all power-suppressed corrections which are required by the Operator Prod- 
uct Expansion (OPE). 

2. The relation of IR renormalons and the OPE has been studied in much detail 
H f|, [5], |(J for the polarization operator of vector currents. In this case the perturbative series 
is complemented by the contribution of the gluon condensate 

Q 2 ^n(Q 2 ) = 1 + ^ + . . . - ^- 2 (9 2 G 2 } + 0(1/Q 6 ) (4) 

By an explicit calculation || it has been shown that the perturbative series diverges, producing 
an ambiguity 0{\/Q 4 ) (s=2,3,. . . in (§)). On the other hand, numerical value of the gluon 
condensate cannot be determined to better accuracy than of order Aq CD because of the quartic 
power divergence. These uncertainties must mutually cancel, since they only arise because of 
our (illegal) attempt to separate perturbative and nonperturbative effects, and an immediate 
question is whether one can organize the expansion in such a way that this problem does not 
appear. A possible solution is suggested by the Wilson OPE, which in fact is not designed to 
separate perturbative and nonperturbative effects, but rather to separate contributions of small 
and large distances. Following this logic literally, we would subtract from the perturbative 
answer the contribution of small virtualities k 2 < fi 2 (where /i is of order 1 GeV), and add it 
to the gluon condensate as a perturbative contribution of order /z 4 . The premium is that the 
separation between the subtracted perturbative answer and the power suppressed contribution 
of the condensate is now unambigous: since the IR region is deleted, the perturbative expansion 
is not plagued by factorially large coefficients and the cancellation of ambiguities between 
perturbative and nonperturbative contributions to the condensate becomes implicit. The price 
to pay, however, is that both perturbative and condensate contributions depend explicitly of 
the scale fi. This dependence may be strong: elimination of the renormalon ambiguity of 
order Aq CD requires reshuffling of a much bigger contribution of order Thus, usefullness 
of this rearreangement should be judged by the practical gain: if renormalon ambiguities 
are numerically much smaller than the phenomenological estimates for the condensates, it is 
hardly reasonable to eliminate the ambiguity at the cost of a large fi dependence. In fact, 
it is easy to see that the whole idea to add HT corrections to the perturbative expansions, 
in which we only know a few first terms and do not see any sign of the factorial divergence 
(and thus expect that we have not yet reached the minimum term), implicitly implies that 
the "true nonperturbative" HT corrections are much bigger than the minimum term in the 
perturbative series, and thus much bigger than the uncertainity in its summation. 

In this respect it is important that semiquantitative estimates of ambiguities in the sum- 
mation of the perturbative series can be worked out. For the particular case of the polarization 
operator using the results of M I get 



a s {Q) 0.002 — 0.02 GeV 4 



l + -^. + ...± 



71 Q 4 



(0.08 - 0.12) GeV 
Q 1 



+ 0(1/Q 6 ) (5) 



where the first number (with an error bar) is an estimate of the renormalon uncertainty^, 

lr The simplest way to make this estimate is by the minimum term in the perturbative expansion. I use a 
somewhat more analytic method and give an imaginary part (divided by tt) of the Borel transform. 



and the second number refers to the phenomeno logical value of the gluon condensate [0]. The 
large range of values for the renormalon ambiguity is due to the fact that it is proportional 
to the fourth power of Aqcd and the uncertainty in the latter is amplified. It is seen that the 
ambiguity is in fact insignificant compared to a 50% error in the phenomeno logical value of 
(g 2 G 2 ). This fact was recognized long ago (see, e.g. ||) and is one of the starting points of the 
QCD sum rules, where the QCD renormalons are ignored as being (supposedly) numerically 
insignificant. 

For the GLS and Bjorken sum rules, assuming the range of values a s (Q 2 = 3 GeV 2 ) = 
0.24 - 0.29, as suggested by the new CCFR data 0, I get 



GLS 
Bj 




a s (Q) ± O02 - 0.07 GeV 2 



7T 



(0.1 ± 0.03) GeV 2 



Ots{Q) 



+ ...± 



0.02 - 0.07 GeV 



7T 



Q 2 



Q 2 

(0.09 ± 0.06) GeV 2 
Q 2 



0(1/Q 4 



(6) 



+ 0(1/Q 4 )(7) 



where the numerical values of the HT 1/Q 2 corrections are taken from the QCD sum rule 
calculations 0. In this case the renormalon ambiguities are roughly factor two smaller than 
the "true" HT effects. Note that the renormalon ambiguities are determined by divergences 
in the coefficient functions in the OPE and are universal in the sense that they do not depend 
on the particular target, while "geniune" HT corrections reflect correlations between partons 
in the target (nucleon) and are process-dependent. 

Finally, it has been shown [|l(| |llj that the pole mass of a heavy quark is affected by IR 
renormalons already at the level of 1/m corrections, which do not allow to determine it from 
the relation to the mass defined at short distances to the accuracy better than 



m pole = m(m) 



3 7T 



± 100 MeV 



This number should be compared to the "true" nonperturbative mass difference between the 
heavy meson and the heavy quark, which is estimated to be of order A = m# — ~ 300 — 500 
MeV 0. 

To summarize, the "renormalon problem" of the separation of perturbative and nonpertur- 
bative contribution can be more or less important in each practical case, depending on which 
method is used to estimate the HT contribution, and what accuracy is claimed. Existing 
calculations of the HT corrections do not pretend to an accuracy better than of order 30-50% 
which is larger (or of order) of the estimated renormalon uncertainty, so that the latter can 
be ignored. However, to improve these predictions one has to treat the renormalon problem 
explicitly. 

3. The IR renormalons have received much attention recently because of their potential 
to expose power-like corrections: the ~ 1/Q 2s uncertainty in summation of the perturbative 
series should be interpreted as indicating presence of the ~ 1/ Q 2s nonperturbative correction 
in the particular physical observable. It should be stressed that this argument cannot compete 
with the OPE if the latter is applicable: attempts to "check" or to "disprove" OPE using the 
renormalon arguments are essentially meaningless. It is precisely the agreement with the 
OPE in its traditional domain which serves as main justification of the hope that the IR 
renormalons may provide a nontrivial information about nonperturbative effects in a more 
general situation, since this technique relies on a purely perturbative analysis and can be 
relatively easily implemented. 



Despite the complicated terminology, the idea of this application is simple, and is that 
IR renormalons in fact probe IR sensitivity of perturbative diagrams to power-like accuracy. 
This can be made explicit using the result of Refs. ||l3t |14| : for IR-safe inclusive quantities 
there is a one-to-one correspondence between IR renormalons and nonanalytic terms in the 
expansion of corresponding Feynman diagrams in powers of some IR regulator like a small 
gluon mass A. In particular, existence of the a/A 2 term in this expansion signals existence 
of the 1/Q uncertainty in summation of the perturbative series, and this necessitates the 
1/Q nonperturbative correction. The A 2 In A 2 term implies the 1/Q 2 correction, etc. For 
illustration, consider the following result for the polarization operator (|4]) calculated with 
a small gluon mass: 



1 + 



o. 



1 



7T 



32 



- 8C(3) 



A^ 

Q 2 



21n(Q 2 /A 2 ) + 



20 



K(3) 



A*_ 

Q 4 



+ 0(A 6 ln 2 A 2 /Q 2 



(9) 



Note that there are no terms ~ A 2 In A 2 (analytic terms like ~ A 2 are not related to the 
IR region) and the first nonanalytic correction is of order ~ A 4 In A 2 , indicating a potential 
nonperturbative contribution of order 1/Q 4 , in agreement with the OPE (|j). Note that absence 
of terms ~ In A 2 is the result of celebrated Bloch-Nordsieck cancellations between contributions 
of real and virtual emission, and absence of certain power-like corrections (alias renormalon 
ambiguities) can be formulated as extension of Bloch-Nordsieck cancellations to power-like 
accuracy [ p~3|j . 

The search of nonperturbative effects using IR renormalons has been been most fruitful in 
heavy quark decays O, O. One finds for the b-quark pole mass 



m 



pole 



m b {m b ) 



4 a s (m) 

3 7T 



A 2 



2 

~a s 

3 nib 



+ 



(10) 



The term ~ indicates presence of an 1/m nonperturbative correction (which is not related 
to matrix element of any local operator and cannot be found using the usual OPE). On the 
other hand, the B-meson total semileptonic inclusive width equals to the same accuracy (for 
simplicity I give the answer for b — > uev transitions, that is for massless quark in the final 
state) 



T(B -> X u ev) 



1927T 3 



m 



pole \ 5 



1-2.41^M + H« s _ 
7r 3 nib 



+ 



(11) 



The ~ vA 2 correction is again present, but is cancelled exactly if the pole mass of the b-quark 
is eliminated in terms of the MS running mass using (|H]). This cancellation presents the 
result of a prime physical importance: inclusive decay widths of heavy particles do not contain 
nonperturbative corrections of order 1/m if they are expressed in terms of mass parameters 
defined at short distances 



10, 13 



Further applications of IR renormalons to the study of nonperturbative efects in resummed 
cross sections will be reviewed by G. Korchemsky ||15|| . 



4. One more idea which has emerged from studies of the QCD renormalons is that 
the Feynman diagrams related to running of the strong coupling, whose low-momentum re- 
gions produce renormalons, may give dominant contributions to perturbative coefficients in 



intermediate orders and can be identified and resummed 0. This can be considered as a nat- 
ural generalization of the proposal by Brodsky-Lepage-Mackenzie (BLM) ]Tf| to eliminate all 
dependence on the QCD /3-function from coefficients of the perturbation theory by adjusting 
the scales of the coupling separately in each order. Several examples show that the second or- 
der perturbative coefficient is significantly reduced by this rearrangement. In high orders one 
again expects a considerable reduction of coefficients since IR renormalons are made implicit 
(they are hidden in uncertainties of the BLM scales |17|). Thus, it is natural to speculate that 
the BLM-improved perturbation theory has smaller coefficients to all orders. 

The corresponding resummation of running coupling effects to all orders has been proposed 
in |17|, [TJ, |14| . Formally, this approach allows to calculate all perturbative corrections of order 

can be traced by contributions of fermion bubble insertions 
into the single gluon line. The resummation of running coupling effects is relatively simple 
and is probably phenomenologically relevant, although the dominance of these correctons is 
a conjecture, which can only be justified a posteriori by comparing to exact calculations. At 
present the corresponding calculations have been done for the heavy quark pole mass [O, WM , 



the Adler function and the r-lepton hadronic width |18|, 14| , and for the exclusive |1| and 



inclusive B-decays. Our results for the r decay show that the value of a s (m T ) extracted 



from these data is probably overestimated, and also give some indication that the commonly 
used resummation of tt 2 corrections is disfavoured in high orders, see |T3[ for details. The 



resummation of (3q corrections in B-decays and its relevance for the extraction of V c b will 
be addressed in her talk by P. Ball |20| . 



5. To summarize, I repeat that the QCD perturbation theory is divergent and does 
not allow to give quantitative predictions to power-like accuracy, unless it is complemented 
by explicit nonperturbative (HT) corrections. In turn, the HT corrections are by themselves 
ill-defined. The corresponding ambiguities have to be in one-to-one correspondence to the 
ambiguities of perturbation theory and must cancel in the sum. For practical cases of the 
GLS and Bj sum rules the corresponding ambiguities are probably a factor 2 below the "true" 
HT corrections. The increasing interest in IR renormalons is trigged by hopes that they can 
help to investigate the structure of nonperturbative corrections in rather general situations, 
and to find physical observables with extended IR stability (to power accuracy), which are 
less sensitive to nonperturbative effects. Another hope is to get estimates for higher-orders 
of perturbation theory, combining the information about the calculated low orders and about 
the renormalons in very high orders. Both directions are interesting and worth efforts. 

As a final remark, let me say that nonperturbative effects in QCD are not reduced to 
renormalons. In particular, not all nonperturbative effects can be traced by divergences of 
the perturbation theory, and also by no means the renormalons can be used to define QCD 
nonperturbatively in the region where the coupling becomes strong. Equally, my talk cannot 
pretend to cover all aspects of QCD renormalons — for example, I ignored ultraviolet renor- 
malons which deserve a special discussion. I refer the readers to the review [^] and original 
papers for the discussion of issues which I was not able to touch here. I thank the organizers of 
this conference for the invitation, and gratefully acknowledge a very rewarding collaboration 
with P. Ball, M. Beneke and V. Zakharov on the subjects related to this talk. 
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